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Abstract — Internet-distributed  hardware-in-the-loop  (ID- 
HIL)  simulation  integrates  geographically-dispersed  HIL 
setups  and  fosters  concurrent  systems  testing  early  in  the 
design  process.  The  degree  to  which  an  ID-HIL  simulation 
loses  fidelity  relative  to  the  single-location  alternative  is 
referred  to  as  distortion.  The  literature  showed  that  various 
factors  affect  distortion,  e.g.,  delay,  jitter,  and  loss.  This 
paper  shows  that  the  coupling  points,  i.e.,  the  particular 
points  at  which  the  system  model  is  divided  to  enable 
distribution  across  the  Internet,  also  affect  distortion.  The 
aim  is  to  turn  coupling  point  selection  into  a  design  decision 
that  can  be  used  to  minimize  distortion.  To  quantify 
distortion,  a  frequency-domain  metric  is  proposed  using  a 
linear  systems  framework.  This  metric  is  then  used  to 
analyze  how  the  choice  of  the  coupling  point  affects 
distortion,  leading  to  guidelines  for  selecting  a  coupling  point 
that  gives  minimal  distortion.  The  theory  is  demonstrated  on 
a  quarter-car  model. 
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I.  INTRODUCTION 

ARDWARE-IN-THE-LOOP  simulation  (HITS)  refers  to 
simulating  a  system  by  coupling  the  physical  models 
of  some  of  its  components  together  with  the  mathematical 
models  of  its  remaining  components  (Bacic  et  al.  2009). 
Thus,  it  combines  the  high  fidelity  of  physical  prototyping 
with  the  cost  effectiveness  of  model-based  simulation 
(Fathy  et  al.  2006).  It  strongly  promotes  concurrent  system 
engineering  and  has  therefore  become  indispensable  in 
many  application  areas,  such  as  automotive  (Kimura  and 
Maeda  1996;  Zhang  and  Alleyne  2005;  Verma  et  al. 

2008) ,  aerospace  (Leitner  2001;  Yue  et  al.  2005;  Cai  et  al. 

2009) ,  manufacturing  (Ganguli  et  al.  2005),  robotics 
(Aghili  and  Piedboeuf  2002;  White  et  al.  2009),  and 
defense  (Huber  Jr  and  Courtney  1997;  Buford  et  al.  2000). 

To  fully  exploit  the  benefits  of  HITS,  it  may  be 
desirable  to  integrate  multiple  HITS  setups  (Keif  2001). 
Recent  efforts  have  focused  on  achieving  such  integration 
over  the  Internet  to  allow  for  integration  of  setups  that  are 
geographically  dispersed  and  unfeasible  to  couple 
physically.  For  example,  the  George  E.  Brown  Jr. 
Network  for  Earthquake  Engineering  Simulation  (NEES) 
(Mahin  et  al.  2003)  provides  an  outstanding  example  of 
the  capabilities  and  impact  of  the  ID-HILS  idea,  and  the 
earthquake  literature  presents  many  other  applications  of 
the  ID-HILS  idea  to  earthquake  simulation  (Watanabe  et 
al.  2001;  Tsai  et  al.  2003;  Spencer  et  al.  2004;  Pan  et  al. 
2005;  Stojadinovic  et  al.  2006;  Mosqueda  et  al.  2008). 
Another  example  in  the  automotive  application  area  is  the 
integration  of  a  ride  motion  simulator  in  Warren,  MI, 
USA,  with  a  hybrid-powertrain-system  simulator  in  Santa 
Clara,  CA,  USA  (Compere  et  al.  2006;  Goodell  et  al. 


2006;  Brudnak  et  al.  2007)  and,  as  a  separate  work,  with 
an  engine-in-the-loop  simulator  in  Ann  Arbor,  MI,  USA 
(Ersal  et  al.  in  press  a;  Ersal  et  al.  in  press  b).  These  works 
highlight  the  potential  impact  of  ID-HILS  for  automotive 
systems. 

Coupling  HIES  setups  over  the  Internet  introduces  a 
deviation  from  the  dynamics  that  would  otherwise  be 
observed  if  the  setups  were  collocated  and  could  be 
directly  integrated.  This  deviation  is  termed  distortion  in 
this  paper. 

There  are  several  sources  of  distortion  in  an  ID-HIL 
setup.  Distributing  a  system  into  subsystems  that  are  co¬ 
simulated  using  independent  numerical  solvers  can  be  in 
and  of  itself  an  important  source  of  distortion  due  to  the 
lack  of  access  to  the  Jacobians  of  the  remote  sites, 
sampling  effects,  etc.,  even  without  any  delay  (Ersal  et  al. 
in  press  b).  Distribution  over  the  Internet  introduces 
further  distortion  due  to  the  Internet’s  delay,  jitter,  and 
loss.  Jitter  refers  to  the  variability  of  delay,  and  loss  means 
that  not  all  packets  sent  arrive  at  destination.  Recognizing 
these  issues,  the  literature  proposed  methods  to  assess  the 
relative  impact  of  distribution  effects  in  comparison  to  the 
effects  of  the  Internet’s  delay,  jitter,  and  loss  (Ersal  et  al. 
in  press  b).  The  literature  also  developed  various 
approaches,  e.g.,  passivity-based  (Anderson  and  Spong 
1989;  Niemeyer  and  Slotine  1991;  Niemeyer  and  Slotine 
2002;  Lee  and  Spong  2006),  event-based  (Xi  and  Tam 
2000;  Elhajj  et  al.  2003;  Mosqueda  et  al.  2008),  and 
observer-based  (Compere  et  al.  2006;  Goodell  et  al.  2006; 
Bmdnak  et  al.  2007)  frameworks,  to  address  stability  and 
distortion  issues  under  a  delayed  coupling  of  subsystems. 

This  paper  focuses  on  another  potential  variable  that 
can  affect  distortion,  namely,  the  coupling  point.  Within 
the  context  of  this  paper,  coupling  point  refers  to  the  point 
at  which  the  HIES  system  is  divided  into  two  subsystems 
that  are  then  co-simulated.  While  options  for  placement  of 
the  coupling  point  may  not  always  exist,  when  they  do 
exit,  the  location  of  the  coupling  point  can  become  a 
design  parameter.  Then,  it  becomes  important  to  know 
how  to  best  pick  that  design  parameter  to  minimize 
distortion. 

Thus,  the  aim  of  this  paper  is  to  develop  a  framework  in 
which  the  impact  of  coupling  point  selection  on  distortion 
can  be  studied  and  the  reasons  that  make  one  coupling 
point  better  than  another  can  be  understood.  As  a  first 
step,  this  paper  will  consider  only  the  delay  and  ignore  the 
jitter,  loss,  and  the  distributed  simulation  effects 
(numerical  issues  due  to  separate  solvers,  sampling 
effects,  etc.).  This  will  not  only  simplify  the  problem,  but 
also  allow  for  the  analysis  to  be  handled  through  a  linear 
framework.  A  linear  framework  will,  in  turn,  allow  for 
leveraging  the  existing  frequency-domain 
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characterizations  of  distortion  (Lawrence  1993;  Yokokohji 
and  Yoshikawa  1994;  ^avu§oglu  et  al.  2002;  De  Gersem 
et  al.  2005;  Griffiths  et  al.  2008). 

The  rest  of  the  paper  is  organized  as  follows.  A 
motivating  example  is  given  first  in  Section  II  that 
illustrates  how  the  location  of  the  eoupling  point  can 
affect  distortion  in  a  quarter-car  representation  of  a 
vehiele.  Then,  in  Seetion  III.A,  a  frequeney-domain 
distortion  metric  from  the  haptics  literature  is  adopted  into 
the  ID-HIL  framework.  Using  this  metrie,  Seetion  III.B 
investigates  which  coupling  point  characteristics  lead  to  a 
low  distortion,  and  relates  distortion  to  a  sensitivity 
function.  Section  III.C  establishes  the  signal-dependence 
of  distortion,  and  Seetion  III.D  discusses  the  effeet  of 
causality  on  distortion.  Finally,  in  Section  IV,  the  theory 
developed  is  applied  to  the  illustrative  quarter-car 
example,  and  eonelusions  are  drawn  in  Section  V. 

II.  Motivating  Example 

Consider  the  system  shown  in  Fig.  1,  which  is  a 
quarter-car  representation  of  a  vehicle,  where  the  masses 
m, ,  nij,  and  represent  the  unsprung,  sprung,  and 
driver  masses,  respeetively,  and  the  eorresponding  spring- 
damper  pairs  capture  the  tire,  suspension,  and  seat 
stiffness  and  damping  properties.  The  figure  also  shows 
the  two  eoupling  point  candidates  eonsidered  in  this  study, 
labeled  as  CPI  and  CP2.  The  coupling  variables  at  these 
eoupling  points  and  their  causality  are  explicitly  shown  in 
Fig.  2.  A  constant  time  delay  of  r  is  considered  in  both 
directions  of  communication,  leading  to  a  round-trip  time 
delay  of  2r  .  The  input  is  the  road  velocity  input,  and  the 
output  of  interest  is  the  displacement  of  the  suspension. 
The  parameter  values  of  the  system  are  given  in  Table  I 
and  are  representative  of  a  military  vehiele  (Ersal  et  al. 
2009). 


Fig.  1  .Example  system  with  two  locations  as  potential  coupling  points 


Table  I 

Parameters  oe  the  Example  System 


Parameter  Value 


*1 

h, 

*3 

K 

ntj 


200  kNs/m 
30  kNs/m 
250  Ns/m 
I  MN/m 
275  kN/m 
1.4kN/m 
1 10  kg 
900  kg 
20  kg 
20  ms 


(a)  (b) 


Fig.  2.  Subsystems  and  coupling  causality  shown  explicitly  for  (a)  CPI 
and  (b)  CP2 

Fig.  3  eompares  the  unit  step  responses  of  the  ideal 
system  and  the  two  systems  in  which  the  coupling 
variables  at  CPI  and  CP2  are  communicated  with  the 
constant  time  delay  of  r  .  As  seen  in  the  figure,  a  delay  at 
CPI  causes  much  more  distortion  than  a  delay  at  CP2. 
This  exemplifies  the  impact  of  the  location  of  the  coupling 
point  on  distortion  and  motivates  the  rest  of  the  paper. 


Fig.  3. Comparing  the  unit  step  responses  of  the  ideal  system  and  the  two 
systems  with  coupling  points  at  CPI  and  CP2 
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III.  Coupling  Point  Analysis 


A.  A  Metric  for  Distortion 

In  this  paper,  distortion  is  a  transfer  function  that  is 
defined  as  the  difference  between  the  reference  dynamics 
Rj  and  actual  dynamics  P ,  where  Rj  represents  the 
dynamics  of  a  HIL  setup  that  is  not  distributed  and  P 
represents  the  dynamics  of  the  same  HIL  setup  distributed 
across  the  Internet.  It  is  also  useful  to  normalize  distortion 
by  the  reference  dynamics  R^ ,  yielding  the  following: 


^  P-Rh 
0  =  -  '' 


R. 


(1) 


This  definition  of  distortion  was  first  introduced  by 
Griffiths  et  al.  within  the  haptics  domain,  where  Rj 
represented  the  dynamics  desired  to  be  rendered  to  the 
user  through  the  haptic  device,  and  P  represented  the 
actual  dynamics  rendered  to  the  user  (Griffiths  et  al. 
2008). 

In  the  following,  an  ID-HIL  system  is  treated  involving 
only  two  sites,  called  local  and  remote.  The  reference 
dynamics  R^  in  this  case  are  achieved  through  an  ideal 
coupling  (involving  bilateral  communications  without 
delay)  of  the  local  and  remote  dynamics.  Fig.  4  depicts  the 
reference  dynamics  in  block  diagram  form,  where  G  and 
G,.  refer  to  the  local  and  remote  dynamics,  respectively, 
u^  is  the  external  input  to  the  local  system,  y,  is  the 
output  of  interest  in  the  local  system,  and  and  are 
the  coupling  variables  between  the  local  and  remote 
systems.  Generally,  the  variables  and  are  power- 
conjugate  variables  modeling  an  energetic  connection, 
such  as  force  and  velocity  in  the  mechanical  domain,  but 
this  need  not  be  true  in  every  ID-HIL  system.  The  desired 
system  equations  are  given  as 


(2) 


Ti 

G,2 

Mj 

T2 

G21 

G22 

U2 

“2  =  Gryi 


from  which  the  desired  dynamics,  R^  ,  from  m,  to  y,  can 
be  derived  as 

„  _G„+(G,G,-G,.G,,)G^  _  G,,+A 


I-G22G, 


I-G,,G, 


where  A=(^G^2G2^-G^lG22)G^ .  Making  the  coupling 
point  explicit  enables  an  analysis  of  the  effects  of 
choosing  different  coupling  points. 


Fig.  4.  Expressing  the  desired  dynamics  in  block  diagram  form 


Next,  to  capture  the  effect  of  delay  due  to  the 
introduction  of  Internet  communications  in  an  ID-HIL 


setup,  consider  a  multiplicative  perturbation.  A,  to  the 
remote  dynamics,  G^ .  This  multiplicative  form  is  suitable 
for  capturing  the  dynamics  of  Internet  delay  and  could 
also  capture  other  unmodeled  dynamics  such  as  the 
dynamics  of  the  sensors  and  actuators.  Fig.  5  expresses  the 
adoption  of  the  distortion  metric  into  the  ID-HIL 
framework  in  block  diagram  form,  where  distortion  is  the 
transfer  function  from  m,  to  . 


Fig.  5.  Adoption  of  the  distortion  metric  into  the  ID-HIL  framework 


The  actual  system  equations  become 


Ti 

Gi2 

_T2. 

G21 

G22 

U2 

=  AG,.y2 


(4) 


where  tildes  are  used  to  differentiate  the  ideal  variables 
from  the  actual  variables.  From  Eq.  (4),  the  actual 
dynamics  P  from  Mj  to  y,  can  be  derived  as 


G„+A^ 

I-AG22G,. 


(5) 


The  distortion  metric  for  the  ID-HIL  framework  can  then 
be  found  as 


A  (1-AG„G,)(G„+/1) 

Eq.  (6)  provides  a  means  to  analyze  the  impact  on 
distortion  of  different  coupling  points  in  an  ID-HIL 
system.  Different  coupling  points  will  lead  to  different 
definitions  of  local  and  remote  dynamics,  i.e.,  different 
G[,,  G[2,  G2, ,  G22,  and  G,, ,  even  though  the  delay 
dynamics  and  other  perturbation  factors  lumped  in  A  may 
remain  invariant,  which  is  assumed  to  be  the  case  in  this 
paper.  Therefore,  different  coupling  points  will,  in 
general,  yield  different  distortion  values,  and  Eq.  (6)  can 
quantify  their  impact  on  distortion. 

B.  Distortion  Analysis 

In  the  framework  created  in  section  III.A,  the  ultimate 
goal  of  bringing  the  actual  dynamics  as  close  to  the 
desired  dynamics  as  possible  translates  to  achieving  a 
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distortion  that  is  as  small  as  possible.  Since  distortion  is  a 
transfer  function  and  thus  a  function  of  frequency,  it  is 
also  possible  to  define  frequency  ranges  over  which  the 
distortion  is  desired  to  be  small. 

Eq.  (6)  reveals  that  there  are  a  number  of  ways  to 
achieve  a  small  distortion  at  a  given  frequency. 
Specifically,  besides  the  trivial  case  of  A  =  1 ,  i.e.,  no 
perturbation,  the  distortion  will  be  small  for  a  given 
frequency,  if  one  of  the  following  is  true  at  that  frequency; 

1.  G,[  — >oo;  the  input  u^  is  greatly  amplified  at  the 
output  y,  and  thus  the  contribution  through  the 
coupling  with  the  remote  system  is  negligible. 

2.  G,2  — >  0  ;  the  feedback  from  the  remote  system 
has  a  very  small  effect  on  the  output  of  interest  y, . 

3.  Gj,  — >0;  the  external  input  u^  has  a  very  small 
effect  on  the  coupling  variable  ■ 

4.  G,.  — >  0  ;  the  remote  system  does  not  affect  the  local 
system,  i.e.,  it  is  driven  by  the  local  system  without 
any  impedance  and  the  coupling  is  almost  one-way. 

5.  G,.  — >oo;  the  remote  system  is  acting  almost  like  a 
wall,  resisting  even  the  smallest  input  . 

6.  G22  — >  ®  ;  the  local  impedance  at  the  coupling  point 
is  so  high  that  the  local  system  is  acting  almost  like  a 
wall  to  the  remote  system. 

Thus,  to  achieve  low  distortion,  one  should  look  for  a 
coupling  point  that  will  lead  to  one  of  the  conditions  listed 
above.  The  physical  interpretation  associated  with  each 
condition  can  help  making  the  assessments  by  inspection. 
However,  it  may  not  always  be  possible  to  observe  that 
one  of  these  conditions  would  be  satisfied  perfectly.  In 
that  case,  it  will  be  necessary  to  consider  all  the  conditions 
simultaneously  when  comparing  coupling  points,  which 
may  make  it  difficult  to  utilize  physical  insight. 
Nevertheless,  although  the  conditions  above  may  seem 
unrelated  at  first  sight,  there  actually  exists  a  concept  that 
unifies  them  all,  and  that  concept  is  the  sensitivity  of  the 
desired  dynamics  to  the  remote  dynamics.  A  formal 
definition  is  given  by 


dGJG^ 


(7) 


In  words,  it  is  defined  as  the  ratio  of  a  relative  change  in 
the  desired  dynamics  to  a  relative  change  in  the  remote 
dynamics.  Evaluation  of  for  the  framework  given  in 
Fig.  4  leads  to 


Sr  = 


GnG.fir 

(1-G22GJ(G,+^) 


(8) 


Comparing  Eq.  (6)  and  (8),  the  following  relationship 
between  distortion  and  sensitivity  to  remote  dynamics  can 
be  derived; 


conditions  and  also  a  single  intuitive,  physical  explanation 
for  distortion.  Also,  this  sensitivity  function  may  be  more 
easily  generalized  to  nonlinear  systems  than  the  distortion 
metric. 

Furthermore,  expanding  the  expression  for  distortion  in 
Eq.  (6)  in  a  Taylor  series  around  A  =  1  shows  that,  to  a 
first  order  approximation,  distortion  is  given  by 
5,  (A -1),  i.e.. 


0  = 


A,(A-l)  +  0((A-lf) 


(10) 


Thus,  to  a  first  order  approximation,  and  recalling  that  A 
is  assumed  to  be  invariant  to  the  location  of  the  coupling 
point,  the  difference  in  distortions  caused  by  different 
coupling  points  is  completely  captured  by  the  sensitivity 
function  .  The  significance  of  this  finding  can  be  seen 
by  referring  to  Eq.  (8)  and  noting  that  can  be  evaluated 
without  the  knowledge  of  the  perturbation  A .  Therefore, 
not  only  provides  a  single  metric  to  be  considered 
when  comparing  coupling  points,  but  also  this  metric  is, 
unlike  0  itself,  is  independent  of  A .  This  allows  for 
comparing  coupling  points  without  having  to  derive  an 
expression  for  A . 

Eq.  (10)  further  implies  that 


^  _SR,IR, 

'■  5G,  / 


(11) 


That  is,  the  sensitivity  is  the  gradient  of  the  distortion 
metric  with  respect  to  the  perturbation  A  at  A  =  1 ,  i.e., 
the  case  when  there  is  no  perturbation. 

Having  related  A,,  to  0  ,  we  can  now  go  back  to  Eq.  (7) 
to  explain  how  a  coupling  point  can  be  selected.  Distortion 
will  be  small  if  a  relative  change  in  the  remote  dynamics 
creates  a  small  relative  change  in  the  desired  dynamics. 
Thus,  the  task  of  finding  the  best  coupling  point  now 
translates  to  finding  the  coupling  point  that  partitions  into 
G,.  all  the  dynamics  whose  relative  change  affects  the 
desired  system  dynamics  the  least. 

C.  On  the  Signal  Dependence  of  Distortion 

It  is  important  to  note  that  distortion  is  defined  for  a 
particular  output  of  interest  y, .  Even  though  the 
formulation  allows  y,  to  be  any  signal  in  the  local  system, 
a  low  distortion  in  y,  does  not  necessarily  imply  that  the 
distortion  will  be  low  in  all  signals  in  the  local  system. 
This  is  easily  demonstrated  by  considering  the  distortions 
in  y^  and  ^2  simultaneously. 

Following  the  same  steps  as  for  y, ,  the  distortion  in  y^ 
can  be  derived  as 


0 


yi 


(12) 


0  =  A, 


I-G22G, 

1-AG„G. 


■(A-1) 


(9) 


From  Eq.  (9)  it  can  be  seen  that  distortion  will  be  small 
when  the  sensitivity  to  remote  dynamics  is  small,  and  it 
can  be  easily  verified  that  the  conditions  listed  previously 
are  also  the  conditions  under  which  becomes  small. 
Hence,  provides  a  unifying  concept  for  those 


Thus,  it  can  be  seen  that,  besides  the  trivial  condition 
A  — >  1  (i.e.,  no  perturbation),  there  is  only  one  condition 
under  which  both  0  and  0^,^  become  small,  namely, 
G,.— >0.  The  conditions  1,  2,  3,  5,  and  6  do  not 
necessarily  imply  a  small  0^^  ,  and  the  condition 
G22  — >  0  ,  which  makes  small,  does  not  necessarily 
make  0  small.  This  emphasizes  the  fact  that  distortion  is 
not  an  independent  property  of  the  system,  but  is  output- 
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signal  dependent,  as  observed  experimentally  before 
(Ersal  et  al.  in  press  b).  Therefore,  when  considering  a 
distortion  analysis,  it  is  important  to  keep  the  signals  in 
mind  with  respect  to  which  distortion  was  defined. 

D.  On  the  Ejfect  of  Causality  on  Distortion 

Under  some  conditions,  distortion  can  be  improved 
simply  by  changing  the  causality.  Consider  the  case  when 
G22  0 .  Furthermore,  assume  that  G22  and  Gf  are 

proper.  In  this  case,  a  switch  in  causality  as  shown  in  Fig. 
6  leads  to  the  following  ideal  system  equations 


metric  for  the  two  coupling  points.  As  the  figure  clearly 
shows,  the  distortion  for  the  system  with  CP2  is  much  less 
than  the  distortion  for  the  system  with  CPI  at  all 
frequencies,  and  hence,  CP2  is  a  better  choice  of  coupling 
point  than  CPI.  The  fact  that  the  distortion  is  on  the  order 
of  1  around  50  rad/s  for  CPI  is  the  dominant  reason  for 
the  significant  difference  between  CPI  and  ideal  response 
in  Fig.  3,  whereas  the  close  agreement  between  the  ideal 
system  and  the  system  with  CP2  is  due  to  the  fact  that  the 
magnitude  of  distortion  is  below  -40  dB  at  all  frequencies 
for  CP2. 


Wj 

U2 

G22 

T2  = 

which  gives 

^  Gi2G2,  +  (1  -  G22G,. ) 

'  G22(1-G22GJ 


(13) 


(14) 


where  the  asterisk  is  used  to  denote  the  switched-causality 
case.  Assuming  a  multiplicative  perturbation  A  as  before 
leads  to 


p-  _  +  G„G22  (a  -  G22G,. ) 

GiilA-GjiGj 


(15) 


Fig.  7.  The  magnitude  of  the  distortion  metric  for  the  two  coupling  points 
CPI  and  CP2 


Finally,  the  expression  for  distortion  is  obtained  as 


^ _ (1-A)G,2G2,G22G, 

(A  -  G22GJ(G,2G2i  +  G„G22  (1  -  G22G,.  )) 

From  Eq.  (16)  it  can  be  seen  that 

lim0*=O  (17) 

G,,-»0 

Thus,  distortion  is  reduced  by  switching  the  causality  at 
the  coupling  point.  Note  that,  in  practice,  such  switching 
may  require  filters  to  ensure  proper  transfer  functions. 


Fig.  6.  Switching  causality  at  the  coupling  point 

IV.  Application  to  the  Example 

In  this  section,  the  theory  that  is  presented  in  Section  III 
is  applied  to  the  quarter-car  example  introduced  in  Section 
II  for  illustration  purposes. 

An  analysis  of  distortion  using  the  proposed  framework 
provides  a  frequency  domain  explanation  into  the  different 
performances  observed  with  the  two  coupling  points  for 
the  same  delay  conditions.  Fig.  7  compares  the  distortion 


To  gain  more  physical  insight.  Fig.  8  compares  the 
sensitivity,  5,. ,  of  the  desired  dynamics  to  the  remote 
dynamics  for  the  two  coupling  points  CPI  and  CP2.  The 
sensitivity  is  less  for  CP2  than  CPI  for  all  frequencies, 
explaining  why  the  distortion  with  CP2  is  lower  relative  to 
the  distortion  with  CP  1 . 

Fig.  9  further  reveals  more  in  detail  why  the  sensitivity 
5,, ,  and  thus  distortion,  is  low  for  CP2  by  showing  the 
magnitudes  of  the  transfer  functions  Gy  ,  i,  j  =  1,2  ,  and 
G,, .  Specifically,  the  distortion  is  low  due  to  Gj2  that 
remains  small  across  the  entire  frequency  range. 
Physically,  this  corresponds  to  the  fact  that  the  driver  and 
seat  subsystem  has  a  very  small  effect  on  the  rest  of  the 
vehicle,  because  the  reaction  forces  of  the  driver  and  seat 
subsystem  have  only  a  small  effect  on  the  sprung  mass, 
which  is  the  largest  mass  in  the  system.  When  this  mass  is 
considered  as  part  of  the  remote  system,  as  is  the  case  with 
CPI,  the  sensitivity  of  the  desired  dynamics  to  the  remote 
system  increases,  which  explains  the  higher  sensitivity  and 
distortion  observed  for  CPI. 

Finally,  Fig.  10  shows  the  magnitude  plot  of  G22  for 
CPI.  The  magnitude  being  significantly  greater  than  zero 
implies  that  distortion  cannot  be  improved  by  reversing 
the  causality  at  CPI.  Thus,  improving  distortion  for  CPI 
requires  other  methods  such  as  feedback  and/or 
feedforward  control.  Development  of  such  methods, 
however,  is  beyond  the  scope  of  this  paper. 
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Fig.  8.  The  magnitude  of  sensitivity  of  the  desired  dynamics  to  the  remote 
dynamics  for  the  two  coupling  points  CPI  and  CP2 


Fig.  9.  The  magnitude  of  G.j,  i,J  =  1,2,  and  G,.  for  CP2  as  factors 
affecting  sensitivity  S,. 


V.  Conclusion 

The  original  contributions  of  this  paper  can  be 
summarized  as  follows.  This  paper  considers  the  coupling 
point  as  a  design  variable  in  ID-HILS  systems  and 
proposes  a  framework  for  a  frequency-domain  analysis  of 
the  impact  of  the  coupling  point  location  to  distortion.  A 
distortion  metric  from  the  haptics  literature  is  adopted  into 
this  framework.  Using  this  framework  and  metric,  the 
paper  identifies  the  system  characteristics  that  render  a 
coupling  point  location  suitable  for  ID-HILS  in  the  sense 
that  it  leads  to  a  low  distortion.  The  paper  further 
identifies  the  sensitivity  of  the  desired  dynamics  to  the 
dynamics  of  the  remote  system  as  the  unifying  reason  for 


different  distortion  results  obtained  with  different  coupling 
points.  It  also  shows  that  distortion  is  an  output-signal 
dependent  concept  and  can,  in  some  cases,  be  affected  not 
only  by  the  location  of  the  coupling  point,  but  also  by  the 
coupling  causality.  The  paper  applies  this  theory  to  a 
quarter-car  model  to  illustrate  and  explain  the  differences 
between  two  coupling  points. 

Future  work  will  focus  on  how  distortion  can  be 
improved  using  feedback  and  how  the  choice  of  coupling 
point  affects  the  trade-off  between  distortion  attenuation 
and  closed-loop  sensitivity.  Other  potential  future 
directions  include  extending  these  results  to  nonlinear 
systems,  as  well  as  considering  the  effects  of  distributed 
simulation,  jitter,  and  loss,  and  the  coupling  of  more  than 
two  subsystems. 
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